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Abstract 



CO 

> 

£C) • When anticommuting Grassmann variables are introduced into a fluid dynam- 



ical model with irrotational velocity and no vorticity, the velocity acquires a 



nonvanishing curl and the resultant vorticity is described by Gaussian poten- 
tials formed from the Grassmann variables. Upon adding a further specific 



interaction with the Grassmann degrees of freedom, the model becomes su- 
persymmetric. 



a.: 

1 I. INTRODUCTION 

> : 

^ ■ An isentropic fluid is described by a matter density field p and a velocity field v. These 

satisfy the continuity equation, which involves the current j = vp, 

p + V-(pv) = (1.1) 

and the force equation 

v + v-Vv = --V? (1.2) 
P 

where P is the pressure. (Over-dot denotes differentiation with respect to time.) We 
show that it is possible to supplement the (p, v) bosonic/commuting variables with Grass- 
mann/anticommuting variables ip such that the entire system exhibits a centrally extended 
super symmetry. Moreover, when the bosonic system is irrotational, so that its vorticity van- 
ishes, Uij = diV-i — djV 1 = 0, and the velocity is the gradient of a velocity potential v = V#, 
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the Grassmann variables give rise to nonvanishing vorticity and provide the Gaussian po- 
tentials in a Clebsch representation for the total velocity (see below). 

The specific system that we analyze devolves from the dynamics for a membrane (a 2- 
dimensional extended object), which propagates in (3 + 1) dimensional space-time. The 
emergent fluid propagates in two spatial dimensions. When the membrane involves just 
bosonic variables, the fluid is irrotational [0]|§. Our supersymmetric fluid is derived by a 
similar construction, starting from a supermembrane [^]. 

In the remainder of this section, we review the action/Hamiltonian formulation for the 
system (|1.1|) - (|1.2|) . In the next section, we directly present the supersymmetric model. 



Section III is devoted to a derivation of this supersymmetric fluid from a supermembrane, 
while concluding remarks comprise the last Section IV. 

For isentropic fluids, the pressure P is a function only of the density, and the right side 
of ( |1.2[ ) may also be written as — VV'(/o), where V is the enthalpy, V"(p) = -P'(p), and 
\f~P' is the sound speed (dash denotes differentiation with respect to argument). Moreover, 
equations (|1 . 1|) and (|1 . 2|) can be obtained by (Poisson) bracketing with the Hamiltonian 



H = j dr^pv 2 + V(p)) (1.3) 
p = {H,p} (1.4a) 
v = {#,v} (1.4b) 
provided the nonvanishing brackets of the fundamental variables (p, v) are taken to be 

K(r),p(r / )} = ^(r-r / ) (1.5a) 

{vHr), v j (r')} = -^^(r - r') . (1.5b) 
P( r ) 

(The fields in the brackets are at equal times, hence the time argument is suppressed.) An 
equivalent, more transparent version of the algebra (|1.5|) is satisfied by the field momentum 
density, which also coincides with the current j. 

V =pv (1.6) 

As a consequence of (|1.5|) we have 

{P i (r),p(r')} = p(r)^(r-r / ) (1.7a) 
{V i (r),V j (r')} = V j (r)d, l 5(r-r')+V i (r')d J 5(r-r') . (1.7b) 
This is the familiar algebra of momentum densities. The Jacobi identity is satisfied by (|1.5| ) 



and ( |1.7| ). The above holds in any dimension |3Q. 

One naturally asks whether there is a canonical 1-form that leads to the symplectic 
structure ( |1.5|) , ( |1.7|) ; that is, one seeks a Lagrangian whose canonical variables can be used 



to derive (|1.5| ) and ( |1 . 7] ) from canonical brackets. When the velocity is irrotational, the 



vorticity vanishes, v can be written as V#, and ( |1.5D is satisfied by postulating that 
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{£(r),p(r')} = 5(r-r') 



(1.8) 



that is, the velocity potential is conjugate to the density, so that the Lagrangian can be 
taken as 



irrotational 



JdrOp-H (1.9) 



where H is given by ( |1.3j ) with v = V0. 

With nonvanishing vorticity, the canonical formulation is more indirect. One writes the 
velocity in a Clebsch decomposition, which in two and three spatial dimensions reads 

v = V6 + aV(3 . (1.10) 

Then 

' drp(6 + a{3) - H . (1.11) 



Here a and (3 are the "Gauss potentials" , and from ( l.llj) is seen that {8, p} as well as 



{P, pa} are canonically conjugate. It then follows that v, given by ( |1.10|) , and p satisfy 

1 Some more observations on the Clebsch decomposition of the vector field v: In three dimensions, 
Eq. ( |1.1C| ) involves the same number of functions on the left and right sides of the equality: three. 
Nevertheless the Gauss potentials are not uniquely determined by v. The following is the reason 
why a canonical formulation of ( |1.5|) requires using the Clebsch decomposition ( |1.10| ). Although 
the algebra (1.5) is consistent in that the Jacobi identity is satisfied, it is degenerate in that the 



kinematic helicity h 



h = \ / d 3 r v • (V x v) = \ / d 3 rv • u> 



(ui l = ^e^ k u>jf.) has vanishing bracket with p and v. (Note that h is just the Abelian Chern- 
Simons term of v.) Consequently, a canonical formulation requires eliminating the kernel of the 
algebra, that is, neutralizing h. This is achieved by the Clebsch decomposition: v = V# + aV/3, 
u = VaxV/?, v-u = V#-(VaxV/3) = V-(#Vax V/3). Thus in the Clebsch parameterization the 
helicity is given by a surface integral h = | / dS • #(Va x V/3) — it possesses no bulk contribution, 



and the obstruction to a canonical realization of (1.5) is removed 



In two spatial dimensions, the Clebsch parameterization is redundant, involving three functions 
to express the two velocity components. Moreover, the kernel of (|ll^) in two dimensions comprises 
an infinite number of quantities 



for which the Clebsch parameterization offers no simplification. (Here u> is the two-dimensional 
vorticity u^j = e^u;.) Nevertheless, a canonical formulation in two dimensions also uses Clebsch 
variables to obtain an even-dimensional phase space. 
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II. SUPERSYMMETRIC FLUID MECHANICS 



A. The Model 

The bosonic fluid model in two spatial dimensions, which descends from a bosonic Nambu- 
Goto action, is supplemented by Grassmann variables ip a that are Majorana spinors (real, 
two- component: ipa = a — 1)2). The Lagrange density reads 



A v 2A 

£ = -p{8 - tyip) - \p{V8 - i^V^) 2 ^— ■ Vip . (2.1) 

P 2 

Here a* are two (i = 1,2), 2x2, real symmetric Dirac "alpha" matrices; in terms of Pauli 
matrices we can take a 1 = a 1 , a 2 = a 3 . Note that the matrices satisfy the following relations, 
which are needed to verify subsequent formulas 

tabalc = e l3 a J ac 
<b°L = 513 5 *c - e lj e ac 

a% ab a% cd = $achd - Sabred + <$aAc (2-2) 



e a b is the 2x2 antisymmetric matrix e = ia 2 . In ( |2.1| ) A is a coupling strength, taken positive. 
The density-dependent potential V(p) = X/p corresponds to a negative pressure P = —2X/p 
and to sound velocity V2X/ p. These describe the "Chaplygin gas". The Grassmann term 
enters with coupling v^2A, so chosen to ensure supersymmetry (see below). It is evident that 
the velocity should be defined as 

v = V6 - l^pVip . (2.3) 

The Grassmann variables directly give rise to a Clebsch formula for v, and provide the Gauss 
potentials. The two-dimensional vorticity reads u = e^d^ = —^e^ditpdjip = — |V?/> x V^>- 



The variables {8, p] remain a canonical pair, while the canonical 1-form in (|2.1|) indicates 
that the canonically independent Grassmann variables are ^/pip so that the antibracket of 
the is 

{^(r),^(r')} = -^(r-r') . (2.4) 



One verifies that the algebra ( |1.5| ) or ( L6|) is satisfied, and further, one has 



{0(r), V(r)} = -— ^(r)5(r - r') (2.5a) 
2p(r) 

{v(r),^(r')} = -^^5(r-r') (2.5b) 

{V(r), ip(r')} = - W(r)5(r - r') . (2.5c) 
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The equations of motion read 

p + V ■ (pv) = (2.6a) 



e + w -ve = \v 2 + ^ + ^^ixx-v^ (2.6b) 



ip + v ■ Vip = ol ■ Vtjj (2.6c) 



and together with (|2.3| ) they imply 

A 

v + v • Vv = V— H Vipct-Vi) . (2.6d) 

P 2 P 

All these equations may be obtained by bracketing with the Hamiltonian 

H = J d 2 r (±pw 2 + + ^ • W) (2.7) 

when (|L5f) , ([Up and (|2.5|) are used. 

We record the components of the energy-momentum tensor, and the continuity equations 
they satisfy. The energy density £ = T°°, given by 



A v 2A 

T°° = ±p W 2 + - + -1— • (2.8) 
P * 

satisfies a continuity equation with the energy flux T°K 

f°° + djT oj = (2.9a) 



T oj = \pv 2 v j - — + ^— Va J v • VV> - -ipd.ip + -e jk iped k ip (2.9b) 
p 2 p p 

This ensures that the total energy, that is, the Hamiltonian, is time-independent. Conserva- 
tion of the total momentum 



d 2 rV (2.10) 
follows from the continuity equation satisfied by the momentum density V 1 = T l ° 

f io + d j T ij = (2.11a) 



9 \ /') X /') \ 

T ij = pv i v j - S ij (— + - — ipa ■ Vip) + - — tpa j diip (2.11b) 



v p 

but the momentum flux T*- 7 , that is, the stress tensor, is not symmetric in its spatial indices, 
owing to the presence of spin in the problem. However, rotational symmetry makes it possible 
to effect an "improvement" , which modifies the momentum density by a total derivative term, 
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leaving the integrated total momentum unchanged (provided surface terms can be ignored) 
and rendering the stress tensor symmetric. The improved quantities are 



(2.12) 



f}° + djTf = 



(2.13a) 



pv % v 3 



P 



-\d k 



(e ki v j + e^v^pipei/j . 
It immediately follows that the angular momentum 

M = I d-r, ij r ; Vj = [ d 2 rpe ij r i v j + \ / dr pc, c 



(2.13b) 



(2.14) 

is conserved. The first term is clearly the orbital part (which still receives a Grassmann 
contribution through v), whereas the second, coming from the improvement, is the spin 
part. Indeed, since |e = \a 2 = S, we recognize this as the spin matrix in (2+1) dimen- 
sions. The extra term in the improved momentum density ^e iJ dj(pipe , ip) can then be readily 
interpreted as an additional localized momentum density, generated by the nonhomogeneity 
of the spin density. This is analogous to the magnetostatics formula giving the localized 
current density j m in a magnet in terms of its magnetization m: j m = V x m. All in all, we 
are describing a fluid with spin. 
Also the total number 



N 



J d 2 r p 



(2.15) 



is conserved by virtue of the continuity equation (|2.6a| ) satisfied by p. Finally, the theory is 
Galileo invariant, as is seen from the conservation of the Galileo boost, 



B = tP - 



J d 2 rrp 



(2.16) 



which follows from (|2.6a|) and (|2.10|) . The generators H, P,M, B and close on the (ex- 
tended) Galileo group. [The theory is not Lorentz invariant in (2 + l)-dimensional space-time, 
hence the energy flux T m does not coincide with the momentum density, improved or not.] 
We observe that p can be eliminated from fl2.1|) so that £ involves only 9 and if). From 
( |2.6b ) and ( |2.6c| ) it follows that 



Substituting into ( |2.1| ) leaves 



£=- 



2A {sj29 -ipip + (V0 - |V^V^) 2 + \il)a ■ Vip) . 



(2.17) 



(2.18) 



Note that the coupling strength has disappeared from the dynamical equations, remaining 
only as a normalization factor for the Lagrangian. Consequently the above elimination of p 
cannot be carried out in the free case, A = 0. 
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B. Supersymmetry 



The theory also possesses supersymmetry. This can be established, first of all, by verifying 
that the following two-component dynamics-dependent supercharges are time-independent 
Grassmann quantities. 



Q a = J d 2 r [pv • (a ab ^ b ) + V2X^a] • (2.19) 



Taking a time derivative and using the evolution equations ( |2.6| ) establishes that Q a = 0. 

Next, the transformation rule for the dynamical variables is found by considering the 
Grassmann charge contracted with a constant Grassmann parameter r] a , giving a bosonic 
symmetry generator Q = t] a Q a . Using the canonical brackets one verifies the field transfor- 
mation rules 

5p = {Q, p} = -V • pivaip) (2.20a) 



\/2A 

56 = {Q, 9} = -±{r)<xip) ■ V9 - foaif)) ■ ipVip + ^— # (2.20b) 

Zp 

\/2A 

5iP = {Q ; ^} = -{rjcx^) ■ Vip - v • ar) rj (2.20c) 



\/2A 

5v = {Q )V } = -(rietip) ■ Vv+ - — rjVifj . (2.20d) 

P 

Supersymmetry is reestablished by determining the variation of the action / dt L, consequent 
to the above field variations: the action is invariant. One then reconstructs the supercharges 
( |2.ig| ) by Noether's theorem. Finally, upon computing the bracket of two supercharges, one 
finds 

{viQa,vlQb} = ^mv2)H (2.21) 

which again confirms that the charges are time-independent: 

{H,Q a } = 0. (2.22) 

Additionally a further, kinematical, supersymmetry can be identified. According to the 
equations of motion the following two supercharges are also time-independent: 



Q a = J d 2 rM. • (2-23) 
Q = ffQ a effects a shift of the Grassmann field: 

(2.24a) 

(#) (2.24b) 

(2.24c) 
(2.24d) 



~6p = 


{Q,p} 


= 


56 = 


{Q,o} 


= -\ 


dip = 


{QM 


= -v 


8v = 


{QM 


= . 
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This transformation leaves the Lagrangian invariant, and Noether's theorem reproduces 
( j2.23j ). The algebra of these charges closes on the total number N. 

{v a iQa,V b 2 Q b } = (viV2)N (2.25) 

while the algebra with the generators fl2.19| ), closes on the total momentum, together with 
a central extension, proportional to volume of space Q — j d 2 r 

{fj a Q a , V b Qb} = {fjart) -P + V2X (fjen)n . (2.26) 

The supercharges Q a ,Q a , together with the Galileo generators (H, P, M, and B), with N 
form a superextended Galileo algebra. The additional, nonvanishing brackets are 

{M,Q a } = \e ab Q b (2.27) 
{M,Q a } = \e ab Q b (2.28) 
{B, Q a } = ot. ab Q b . (2.29) 



III. MEMBRANE CONNECTION 

The equations for a supersymmetric Chaplygin fluid devolve from the supermembrane 
Lagrangian, Lm- We shall give two different derivations of this result, which make use of 
two different parameterizations for the parameterization-invariant membrane action and give 



rise, respectively, to fl2.1|) and (|2T 



We work in a light-cone gauge-fixed theory: The membrane in 4-dimensional space- 
time is described by coordinates x^ (fi = 0,1,2,3), which are decomposed into light-cone 
components x = ^(x° ±x 3 ) and transverse components x % {% = 1, 2}. These depend on an 
evolution parameter r and two space-like parameters <p r {r = 1,2}. Additionally there are 
two- component, real Grassmann spinors ip, which also depend on r and cff . In the light-cone 
gauge, x + is identified with r, x~ is renamed 8, and the supermembrane Lagrangian is 

L M = [d 2 (f>£ M = - [ d 2 (f>{VG-le rs d r iPcxd s 4)-x} (3.1) 

where G = det G a /3', 



G 




(3.2) 



a/3 — 

2d T 9 — <9 r x ■ <9 r x — i/jd T vjj u s 

U r Q rs 

G = gT 

T = 2d T 9 — 9 T x • <9 r x — ipdrip + g rs u r u s 
g rs = 9 r x ■ <9sX , g = det g rs 

u r = d r 6 — \il)d r ip — <9 T x ■ 9 r x . (3.3) 



Here d T signifies differentiation with respect to the evolution parameter r, while d r ,d s dif- 
ferentiate with respect to the space-like parameters (<f) r , <p s ) , and g rs , the inverse of g rs , is 
used to move the (r, s) indices. Note that the dimensionality of the transverse coordinates 
x l is the same as of the parameters r , namely two. 



A. First Derivation 



To give our first derivation, we rewrite the Lagrangian in canonical, first-order form, with 
the help of canonical momenta defined by 

tttt - P — n<9 r x — Hu r <9 r x (3.4a) 
oo T x 

dC M 



C M = p • c> T x + Ud T 9 - ^d T ij + 7 ^(p 2 + g) + ±e rs a r ^£*<9^ ■ x 

+ u r (d r x ■ p + Ud r 6 - \X\ij)d r if) . (3.5) 

In ( |3.5| ) if serves as a Lagrange multiplier enforcing a subsidiary condition on the canonical 
variables. The equations that follow from ( |3.5|) coincide with the Euler-Lagrange equations 
for (|3.1| )- (|3.3j ). The theory still possesses an invariance against redefining the spatial pa- 
rameters with a r-dependent function of the parameters. This freedom may be used to set 
u r to zero and fix II at —1. Next we introduce the hodographic transformation [||], whereby 
independent-dependent variables are interchanged, namely we view the cff to be functions 
of x l . It then follows that the constraint on (|3.5|), which with II = —1 reads 



<9 r x ■ p — d r 9 + ^ipd r ip = (3.6a) 

becomes 

«9 r x . ( p - V0 + §VVV>) = . (3.6b) 

Here p, 9 and ip are viewed as functions of x, renamed r, with respect to which acts the 
gradient V. Also we rename p as v, which according to ( |3.6b| ) is 



v = V6 - . (3.6c) 

From the chain rule, it follows that 

d T = d t + <9 r x ■ V (3.7) 

and according to ( |3.4a| ) (at n = — 1, u r = 0) <9 T x = p = v. Finally, the measure transforms 
according to d 2 (f) —> d 2 r -^=. Thus the Lagrangian for ( |3.5| ) becomes, after setting if to zero 
and n to —1, 
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J y 9 

But e TS dgxi d r x % = e u det <9 r x J = e^-^/g. After is renamed \/2A/p, fl3~8|) finally reads 

1 /■ A \ /2 A 
7 =) y d 2 r - - §p(V0 - \ipVif)) 2 x VV) • (3.9) 



Upon replacing ^ by ^(1 — e)^, this is seen to reproduce the Lagrange density ( |2.1| ), apart 
from an overall factor. 



B. Second Derivation 



For our second derivation, we return to ( p. 1| ) — ( pT3|) and use the remaining reparameter- 
ization freedom to equate the two x l variables with the two r variables, renaming both as 



r % ||. Also t is renamed as t. In (p. 1|) ( pT3|) g TS = 5 rs , and <9 T x = 0, so that ( |3.3| ) becomes 
simply 

G = T = 26 -ifjij) + u 2 (3.10) 

n=V6 -\if)Vip . (3.11) 



Therefore the Nambu-Goto action (3.1) reads 



L M = - J d\y28 - ifnp + (Ve - \1pV1P) 2 + \ificx x V^} . (3.12) 

Again a replacement of if) by -^(1 — e)if> demonstrates that the integrand coincides with the 
Lagrange density in ( |2.18| ) (apart from a normalization factor). 

C. Further Consequences of the Supermembrane Connection 

The supermembrane dynamics is Poincare invariant in (3+l)-dimensional space-time. 
This invariance is hidden by the choice of light-cone parameterization: only the light-cone 
subgroup of the Poincare group is left as a manifest invariance. This is just the (2 + 1) 
Galileo group generated by if, P, M, B, and N. (The light-cone subgroup of the Poincare 
group is isomorphic to the Galileo group in one lower dimension.) The Poincare generators 
not included in the above list correspond to Lorentz transformations in the "— " direction. 
We expect therefore that these generators are "dynamical", that is, hidden and unexpected 
conserved quantities of our supersymmetric Chaplygin gas, similar to the situation with the 
purely bosonic model 0. 

One verifies that the following quantities 

D = tH-Jd 2 rp6 (3.13) 

G = y d 2 r (rS - VP i - \if)aot ■ Vj$) (3.14) 
= Jd 2 r(r£ - 6P - \if)cx.ot ■ Pip) (3.15) 
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are time-independent by virtue of the equations of motion (|2.6| ), and they supplement the 
Galileo generators to form the full (3+1) Poincare algebra, which becomes the super-Poincare 
algebra once the supersymmetry is taken into account. 

IV. CONCLUSION 

We have shown how fluid dynamics can be extended to include Grassmann variables, 
which also enter in a supersymmetry-preserving interaction. Since our construction is based 
on a supermembrane in (3+l)-dimensional space-time, the fluid model is necessarily a planar 
Chaplygin gas. It remains to be shown how this construction could be extended to arbitrary 
dimensions and to different interactions. Note that Grassmann Gauss potentials can be used 
even in the absence of supersymmetry. For example, our theory fl2.1p , with the last term 
omitted, posseses a conventional, bosonic Hamiltonian without supersymmetry, while the 
Grassmann variables are hidden in v and occur only in the canonical 1-form. In a related 
investigation, conventional fluid mechanics is generalized, so that it possesses a non-Abelian 
gauge symmetry ||. 

Note Added: J. Hoppe has informed us that some of the above results were obtained 
by him in unpublished research: Karlsruhe preprints KA-THEP-6-93 and KA-THEP-9-93 
( |hep-th/931i059D . 
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